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ABSTRACT 

We investigate the effect of blending in future gravitational microlensing surveys by carrying out simulation 
of Galactic bulge microlensing events to be detected from a proposed space-based lensing survey. From this 
simulation, we find that the contribution of the flux from background stars to the total blended flux will be 
equivalent to that from the lens itself despite the greatly improved resolution from space observations, implying 
that characterizing lenses from the analysis of the blended flux would not be easy. As a method to isolate 
events for which most of the blended flux is attributable to the lens, we propose to use astrometric information 
of source star image centroid motion. For the sample of events obtained by imposing a criterion that the 
centroid shift should be less than three times of the astrometric uncertainty among the events for which blending 
is noticed with blended light fractions / B > 0.2, we estimate that the contamination of the blended flux by 
background stars will be less than 20% for most (~ 90%) of the sample events. The expected rate of these 
events is > 700 events/yr, which is large enough for the statistical analysis of the lens populations. 
Subject headings: gravitational lensing 



1. INTRODUCTION 

Observational experiments to detect light variations of 
stars induced by gravitational microlensing h ave been an d 



are b eing carri ed out for more than a d ecade ( Udalski et al 
1993| lAlcock et alJ fi99l lAubourg et alJ fi993t iBond et al' 
2001). Originally, the experiments were initiated for the 



purpose of sear ching for dark matter in the Galactic halo 
JPaczvriskHll986l) . However, microlensing developed many 
applications in various aspects of galactic, stellar, and plan- 
etary astrophys ics such as probing the structure of the inner 
gala xy ( Kirag a & PaczvriskJl99^lEvansll994tlHan & Gouldl 
1 1 1995bl lZhao. Spergel. & Rich! 11995ft resolving the pho- 
tosphere"^ dis t ant stars dWitll 119951 lAlbrow et all 119991: 
lAbe et alJ 120031 IField et alJ 120031), a nd detecting e xtraso- 
lar planets~llMao & Paczvris kil 1 1 99 It iGould & Loebl [ 1992 1 : 
[ Bond et alJ |2004t lUdalski et all 120051 iBeaulieu eta!! hood 
Gould et al. 2006). Current experiments are routinely detect- 
ing more than 500 events every season and the total number 
of detected events is now reaching 3000. 

The probability of a star being gravitationally lensed is 
very low (of the order of 10~ 6 ). As events are very rare, 
lensing searches are conducted toward very crowded star 
fields such as the Galactic bulge in order to monitor as many 
stars as possible within a single frame. For observations in 
such a dense field, blending of stellar images is inevitable 
SDi Stefano & Esidll995tlWozniakll 19971) . The lightcurve of 
a blending-affected event is expressed as 



Fobs = AFq + Fb, 



A = 



u 2 + 2 
k(w 2 +4)' 



(1) 



where Fq is the baseline flux of the lensed source star, Fg is 
the amount of the blended flux, A is the magnification induced 
by lensing, and u is the lens-source separation normalized by 
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the Einstein ring radius f?E- The normalized lens-source sepa- 
ration is related to the lensing parameters by 



it = 



t-tp 

tE 



1/2 



(2) 



where fg is the Einstein timescale, fo is the time of the closest 
lens-source approach, and uq is the separation at that moment 
(impact parameter). Among the lensing parameters, the Ein- 
stein timescale is related to the mass and location of the lens 
and relative lens-source transverse speed, and thus it can pro- 
vide constraint about the lens system. However, the constraint 
is weak because the timescale results from the combination of 
the physical parameters. Blending aggravates the photometric 
precision due to the increased noise from the blended flux. In 
addition, blending causes accurate determinations of the lens- 
ing parameters difficult due to the poorly known baseline flux 
of the source star. As a result, the constraint about the lens 
system becomes further weaker. 

Recently, next-generation microlensing experiments based 
in space are being seriously discussed. Microlensing Planet 
Finder (MPF), that succeeds the original con cept of Galactic 
Exoplanet Survey Telescope (GEST) mission (Ben nett & Rhid 
2002), is a space mission proposed to NASA's Discovery Pro- 
gram with the main goal of searching for a large sampl e of ex- 
trasol ar planets by using gravitation lensing method JBennettl 
12004ft . By conducting observations in space, it is possible to 
minimize blending effect thanks to the improved resolution. 
Another important advantage of the space observation may be 
that one can directly detect and characterize the lens based on 
the observed flux. This is possible when a light-emitting star 
is responsible for a lensing event and the flux from the lens 
contributes a substantial fraction of the observed flux. It is 
believed that majority of e vents detected tow ard the Galactic 
bulge are caused by stars (lHan & G ould 2003). In addition, 
the MPF mission plans to monitor faint main-sequence stars 
in order to optimize detections of low-mass planets by mini- 
mizing finite-source effect. Therefore, the light from the lens 
would not be negligible for a significant fraction of the event 
sample, and thus the existence of the lens could be revealed 
by the excess flux from the lens. By conducting multi-band 
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observations, it would be possible to determine the spectral 
type of the lens, which would provide much tighter constraint 
about the lens. 

However, this new method of constraining lenses can be 
applicable under the condition that source stars to be ob- 
served are in most cases free from blending by the flux of 
blends other than lenses. If an unrelated background star is 
located within the resolution disk of the lensed source star, the 
blended flux comes not only from the lens but also from the 
background star. This causes complexity in the analysis of the 
observed blended light , makin g it difficult to uniquely char- 
acterize the lens. iHanl l|2005) (hereafter Paper I) warned of 
the possible contamination of the blended flux by background 
stars. As a method to resolve the contamination, he suggested 
to use the astrometric information of the source star image 
centroid shift. If an event is affected by background stellar 
blending, the centroid will shift from the apparent position of 
the combined image of the sou rce plus blended stars toward 
the lensed star during the event (Goldberg 1998). 

In this paper, we further investigate the blending problem 
in future space-based lensing surveys following Paper I and 
probe the feasibility of sorting out events whose blended flux 
is mostly attributable to the lenses by using the additional 
astrometric information of the centroid shift. For this, we 
carry out simulation of Galactic bulge microlensing events 
expected to be detected from the proposed MPF survey by 
imposing realistic observational conditions and detection cri- 
teria. Based on this simulation, we investigate the expected 
pattern of blending and the seriousness of the contamination 
by background stars. We also estimate the fraction of events 
that can be identified to be free from background stellar blend- 
ing by using the proposed method. 

2. SIMULATION 

The basics scheme of the simulation of lensing events is 
similar to the one in Paper I. The locatio ns of the lenses 
and source stars are allocated based on lHan & Gould fl995a) 
Galactic mass distribution model, that is composed of a 
double-exponential disk and a barred bulge. The velocity dis- 
tribution of the bulge is deduced from the tensor virial theo- 
rem, while the disk velocity distribution is modeled to have 
a flat rotation speed of v c = 220 km s" 1 and velocity disper- 
sion along and normal to the disk plane of cni = 30 km s -1 
and <7j_ = 20 km s" 1 , respectively. The brightnesses of the 
bulge stars are assigned on the basis of a combined luminos- 
ity fun ction of constructed based on t hose o f iHoltzman et al.l 
(1998) and lGould. Bahcall. & Frvnnl (11991) considering ex- 
tinction and distance modulus. The lens masses are as- 
signed based on the mass function of lGouldl (2000), which 
is composed of stars and stellar remnants. The mass function 
model is constructed under the assumption that bulge stars are 
formed initially according to a double power-law distribution 
oidN/dM = k(M/0J MqV, where 7 = -2.0 for M > 0.7 M© 
and 7 = -1 .3 for M < 0.7 M©. Based on this initial mass func- 
tion, remnants are modeled by assuming that the stars with 
initial masses 1 M© < M < 8 M©, 8 M© < M < 40 M©, and 
M > 40 M© have evolved into white dwarfs (with a mean 
mass (M) = 0.6 M©), neutron stars ((M) = 1.35 M©), and 
black holes ((M) = 5 M©), respectively. For events caused 
by stellar lenses, the lens brightness is determined based on 
the lens mass by using a mass-luminosity relation. For events 
caused by remnant lenses, the lens is assumed to be dark. 

However, the main focus of Paper I was investigating pos- 
sible contamination of blended flux in general and thus it is 



not focused on a specific survey. As a result, there are several 
major differences between the simulations in Paper I and this 
work. The differences are described below. 

1 . Since no specific instrument was defined, the analysis 
in Paper I was based on all possible events regardless of 
their detectability. In this work, analysis are based only 
on events that could be detected from the MPF survey. 
The detectability of an events is subject to the condition 
and strategy of observation. For the observation field, 
we assume that the future survey will monitor toward 
a field centered at (l,b) ~ (1.2°, -2 .4°) with a field of 
view of (A/, A/?) - (0.93°, 2.8°). For a sampling fre- 
quency, we assume that the survey will continuously 
observe the target field with an interval of 15 minutes 
during 9 months of a year. For each interval, we as- 
sume that a 10-minute exposure image is acquired with 
a photon count rate of 1 3 photons/sec for an / = 22 star. 

2. Another factor that determines the event detectability is 
the selection criteria. We judge the detectability of an 
event based on the uncertainties of the lensing param- 
eters recovered from the lightcurve. We determine the 
uncertainties by computing the curvature matrix of % 2 
surface of the lensing parameters for each event pro- 
duced by the simulation. For the case of a lensing 
lightcurve, the curvature matrix is defined as 



A dFabs t dF ( 



"obs.k 



1 



dpi dpj a\ 



(3) 



where Af bs is the number of observations, crj = \/F b s 
is the photometric precision of each measurement, and 
p, = (Fo,FB,uo,to,tE) are the five lensing parameters re- 
quired to fit the lightcurve of a standard point-source 
single-lens event. Then, the uncertainties of the indi- 
vidual lensing parameters correspond to the diagonal 
components of the inverse curvature matrix (covariance 
matrix), i.e. 

c = b~ l . (4) 



Cu, c — b . 

With the determined uncertainties, we set the detection 
criteria of events such that the fractional uncertainties 
of the timescale and blended flux should be less than 
20%, i.e. a, E /t E < 0.2 and a FB /F obs < 0.2. Under this 
definition of event detectability, we note that events can 
be detected even when the source star trajectory does 
not enter the Einstein ring of the lens. This can be seen 
in the upper left panel of Figure^ where we present the 
distribution of impact parameters of events. From the 
distribution, one finds that a fraction of events have im- 
pact parameters larger than unity, especially for events 
involved with relatively bright source stars (Iq < 22). 
In Paper I, detectable events were limited only to those 
with source trajectories entering the Einstein ring, i.e. 
uo < 1 .0. 

3. Another major difference is the brightness range of 
source stars. In Paper I, we considered only stars with 
/-band baseline brightness of Iq < 23. With the im- 
plementation of the detection criteria described above, 
we found that events associated with much fainter stars 
could be detected. In this work, we, therefore, do not 
set the limit of source star brightness as long as the 
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TABLE 1 
MPF Event Rate 



categories 



20 22 24 26 

source star brightness, I 

FIG. 1. — Distribution of lensing events in the parameter space of the 
/-band baseline source brightness, Iq, and lens-source impact parameter, uq. 
The distribution is for events to be detectable from a lensing survey equipped 
with a space telescope similar to the Microlensing Space Finder mission. The 
contour levels are in relative scale. 



event meets the detection criteria. With the new cri- 
terion, we find that events can be detected with source 
brightness down to Iq ~ 27, as shown in the lower panel 
of Figure [2 where the distribution of /-band baseline 
source star brightness is presented. With this new range 
of the source star brightness, the baseline flux of the 
source star decreases in average while the blended flux 
remains similar to the previous estimation. As a re- 
sult, the average value of the blended light fraction, 
/b = Fs/(Fq+Fb), is smaller than the value in Paper I. 

The last major change is the definition of the resolution 
angle f? res . In Paper I, we assume that two stars cannot 
be resolved if the separation between them is less than 
a fixed value of the diffraction limit (DL) regardless of 
the flux ratio between the two stars. To describe the ex- 
periment more realistically by considering the variation 
of the resolution angle depending on the relative flux, 
we set the resolution angle as 



0.5 FWHM 

[-0A4(F 2 /Fi) + 0.54] FWHM 



if F 2 /Fi < 0.1, 
otherwise, 



(5) 

where F\ and F2 are the fluxes of the brighter and fainter 
stars, respectively, and FWHM ~2x DL ~ 0.4" corre- 
sponds to the diameter of the resolution disk. Then, 
the resolution angle for two equally bright stars is # res = 
0.1 FWHM - 0.04" 4 With the adoption of the new 
definition of the resolution angle, the average blended 

4 For the MPF mission, the pixel size of the detector is ~ 0.24"/pixel, 
which is bigger than the assumed DL = 0.2". However, we note that images 
to be taken by the MPF will be dithered and there will be many of them. 
Therefore, it would be possible to create drizzled images that recover the 
diffraction limit. 



total event rate 

- bulge self-lensing events 

- disk-bulge events 

- stellar lens events 

- remnant events 



3415/yr 

2192/yr (64.2%) 
1223/yr (35.8%) 
2510/yr (73.5%) 
905/yr (26.5%) 



NOTE. — Detection rate of microlensing events 
from the MPF survey. 



light fraction further decreases compared to that in Pa- 
per I. 

Based on the simulation, we estimate that the total event 
rate of the MPF survey would be r tot ~ 3500 events/yr, 
among which bulge self-lensing and disk-bulge events com- 
prise ~ 64% and ~ 36%, respectively. The ratio between 
the events caused by stellar and remnant lenses is I\ : r rem = 
73.5 : 26.5. We summarize the result of event rate in Tabled 

3. RESULTS 

The top panel of Figure |2] shows the distribution of events 
expected to be detected from the MPF survey in the parameter 
space of the Einstein timescale and blended light fraction. We 
compare the distribution with those of events expected when 
the only blending source is the lens (middle panel) and back- 
ground stars (bottom panel). In each panel, the points marked 
by bright and dark-tone greyscales represent events caused by 
stellar and remnant lenses, respectively. 

In the case where the dominant source of blending is the 
lens itself, the distribution has the following distinctive fea- 
tures. 

1. First, the region with very little blended light fraction 
(/b ^ 0.05) is densely populated by events caused by 
remnant lenses, which are dark and thus have no con- 
tribution to the blended flux. 

2. Second, except for the ones caused by remnants, events 
are distributed smoothly over the entire range of the 
blended light fraction. 

3. Third, there exists a correlation between the timescale 
and blending fraction. This correlation arises because 
heavier lenses tend to be brighter and thus events with 
longer timescales are more likely to be affected by 
larger amount of blended flux. We find that the mean 
timescale is (fg) ~ 20 days for stellar-lens events with 
/b < 0.1, while (r E ) — 35 days for events with / B > 0.9. 
However, this correlation is not very strong due to the 
large dispersions of the Einstein timescale and lens 
brightness. 

In the case where most of the blended flux comes from 
background stars, on the other hand, the distribution has very 
different features from those described above. These features 
are as follows. 

1. First, the distribution is not smooth, but divided into 
two groups of heavily and lightly blended events. The 
two groups are roughly divided by /b ~ 0.5. This kind 
of distribution occurs when the number of blended stars 
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lens + background stars 
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FIG. 2. — Distribution of microlensing events detectable from the MPF sur- 
vey in the parameter space of the Einstein timescale, ?e. and the blended light 
fraction, /b (top panel). The other two panels are the distributions expected 
when the only blending source is the lens (middle panel) and background 
stars (bottom panel), respectively. In each panel, the points marked by bright 
and dark-tone greyscales represent the events caused by stellar and remnant 
lenses, respectively. 



in the resolution disk is of an order of one or less. If the 
number of blended stars is substantially larger (smaller) 
than this, most events will be distributed in the region 
of high (low) blended light fraction fraction. We find 
that the average number of blended stars is ~ 0.42. 

2. Another noticeable feature is that both of the events 
caused by stellar and remnant lenses have similar dis- 
tributions. This is because there is no contribution of 
the flux from the lens to the total blended flux. 

Considering the characteristics of the distributions for the 
stellar and lens blending dominated regimes, we find that the 
distribution of MPF events has features of both types of blend- 
ing. For example, although majority of remnant events are 
located in the very small /b region, a considerable fraction 
of these events are smoothly distributed in a wide range of 
blended light fraction. In addition, although still events can be 
divided into two groups of heavily and lightly blended events, 
the distinction between the two groups is less clear. These 
characteristics imply that both types of blending are impor- 
tant for events to be detected from the MPF survey. This can 
be seen also in Figure [3] where we plot the distribution of 
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FIG. 3. — Distribution of blended light fractions, /b, for microlensing 
events to be detectable from the MPF survey (solid curve). The dotted and 
dashed curves are the distributions when the only blending source is assumed 
to be the lens and background stars, respectively. Lower panel shows the 
cumulative distributions. 



blended light fraction. From the figure, one finds that both 
types of blending are nearly equally important. 

4. ADDITIONAL ASTROMETRIC INFORMATION 

When a source star blended with background stars is grav- 
itationally magnified, the apparent position of the blended 
image centroid moves toward the source position during the 
event. The amount of the centroid shift is 



Ar=|r-ro|=X>|rs-r |; V = 



(1-/ B )(A-1) 
(1-/b)(A-1)+1 ! 



(6) 



where rs is the position of the source, and fo and r are the 
positions of the image centroids before and in the middle of 
lensing magnification, respectively. 

The centroid shi ft can be m easured in other way. iHanl 
(2000) and lGould & Anl (120021) pointed out that the position 
of the lensed source star, rs, instead of the centroid position of 
the combined image of the source and blend, f , can be directly 
measured on an image obtained by using difference imag- 
ing. Difference imaging is an image-subtraction technique to 
provide accurate photom etry and astrometry of variable stars 
in crowded fields (Tomanev & Crotts 1996; Alard & Lupton 
1998). In this technique, one first forms a high-quality "tem- 
plate" image. For each of other images ("current" image), 
one convolves the template image to the same seeing as the 
current image, translates it so that the two images are geomet- 
rically aligned, and linearly rescales its flux so that the two 
images are photometrically aligned as well. When the con- 
volved template image is subtracted from the cunent image, 
all that remain are a set of point spread functions (PSFs) at the 
locations of stars that experienced light variation. Then, by 
transforming the coordinates of the source position measured 
on the subtracted image to the template image, one can mea- 
sure the positional difference between rs and Fq. The amount 
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FIG. 4. — Distribution of source star image centroid shifts, Adi for mi- 
crolensing events to be detectable from the MPF survey. The dotted and solid 
curves are for all events and those with blended light fraction /rj > 0.2, re- 
spectively. Lower panel shows the cumulative distributions. 



of the centroid shift measured in this way is related to the shift 
measured directly on two regular images by 



Ar DI = |r s -r | 



(7) 



We note that the centroid shift measured in this way is always 
larger than the shift in equation © because T> < 1.0. 

Astrometric information of the centroid shift might be im- 
portant in distinguishing events blended by the lens and back- 
ground stars. This is because typical separation between the 
source and blended background star, which is of an order of 
10~ 2 arcsec, is much larger than the separation between the 
source and lens, which is of an order of 10~ 4 arcsec. Then, 
if no centroid shift is detected within the astrometric preci- 
sion of the centroid shift measurement despite that blending 
is noticed from the lightcurve, the chance that the lens is the 
blending source is very high. Figure0]shows the distribution 
of the centroid shifts for events detectable by the MPF survey. 

We investigate the usefulness of the astrometric informa- 
tion in sorting out events where most of the blended flux 
is attributable to the lens. For this, we calculate the frac- 
tion of blended flux from background stars, F M out of the 
total blended flux for events where blending is noticed with 
/b > 0.2 but no centroid shift is detected within the astromet- 
ric uncertainty. The astrometric uncertainty of the centroid 
shift measurement is 



FWHM 



fast ' 



1/2 



(8) 



where N v is the photon count. We assume that the centroid 
shift is measured from difference imaging, i.e. An = |rs — fo|. 
Despite the absence of atmosphere, difference imaging on 
space-based observation may not be simple due to variations 
in the PSF profile from image to image cause d by temporal 
variat ions in the thermal load on spacecraft ( Gilliland et al 
2000). Furthermore, the MPF mission may suffer from un- 



0.5 
0.4 
0.3 
0.2 
0.1 
1 

0.8 
0.6 
0.4 
0.2 




for ^vents with / B ^0.2 
! w/o astrometric information 

i — with centroid shift §3(7 . 




,90% 



4 9'K 



0.2 



0.4 0.6 

F/F a 



FIG. 5. — Distribution of the fractions of the blended flux from background 
star, F*, out of the total blended flux, Fb, for MPF microlensing events with 
the blended light fraction /b > 0.2 (dotted curve). Solid curve is the distribu- 
tion for events that are additionally filtered out with astrometric information 
under the condition that the measured source star image centroid shift, Arm, 
is less than 3 times of the astrometric uncertainty, cr as t ■ Lower panel shows 
the cumulative distributions. 



dersampling of the PSF which further complicates the ge- 
ometric alignment procedure which involve interpolation of 
undersampled PSFs. To account for these possible additional 
sources of uncertainty in astrometric measurement, we con- 
sider only centroid shifts with Aroi > 3(7 ast can be firmly de- 
tected. 

In Figure [5] we present the resulting distributions of the 
fraction F^/Fr. Dotted curve is for events with blended light 
fractions f B > 0.2 and solid curve is the distribution for events 
that are additionally filtered out with astrometric informa- 
tion under the condition that the measured centroid shift is 
ArDi < 3<7 as t. We note that the blended flux from background 
stars works as a contaminant to the flux from the lens, and thus 
small ratio of F+/Fb implies that the blended flux is mostly 
from the lens. From the distributions, we find that blended 
flux of nearly half (49%) of all events with /b > 0.2 is con- 
taminated by the flux from background stars by more than 
F±/Fb = 20%. However, for the sample of events filtered out 
by using the additional astrometric information, it is found 
that only a minor fraction (~ 10%) of events are contaminated 
by background stellar blending. We also find that this sample 
of events comprises ~ 22% of all MPF lensing events. Con- 
sidering the total event rate of ~ 3500 events/yr of the survey, 
the rate of these events is > 700 events/year, which is large 
enough for the statistical analysis of the lens populations. 

5. CONCLUSION 

We investigated the effect of blending in future gravitational 
microlensing surveys using aim class space telescope. For 
this, we carried out simulation of Galactic bulge microlensing 
events to be detected from the proposed MPF lensing survey 
by imposing realistic observational conditions and detection 
criteria. From this simulation, we found that the contribu- 
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tion of the flux from blended background stars to the total 
blended flux is equivalent to that of the lens itself, implying 
that characterizing lenses from the analysis of the blended flux 
would not be easy. However, with the additional astrometric 
information of the source star image centroid shift, we found 
that it would be possible to isolate events for which most of 
the blended flux is attributable to the lens. For the sample 
of events obtained by imposing a criterion that the centroid 
shift should be less than three times of the astrometric uncer- 
tainty among the events for which blending was noticed with 



blended light fractions /b > 0.2, we estimated that the con- 
tamination of the blended flux by background stars would be 
less than 20% for most (<~ 90%) of the sample events. The ex- 
pected rate of these events is > 700 events/yr, which is large 
enough for the statistical analysis of the lens populations. 
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